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The  Disk:  A  Comparison  of  Electromagnetic 
Scattering  Solutions  and  Its  Use  as  a 
Calibration  Standard  for  Bistatic  RCS  Measurements 


1.  INTRODUCTION 


This  report  compares  three  circular  disk  solutions:  a  closed  form  physical 
optics  (PO)  solution,  the  RADC/Syracuse  Research  Corp.  bistatic  radar  cross- 
section  (BSRCRCS)  computet  program  that  utilizes  the  physical  theory  of  diffraction 
(PTD),  and  an  exact  eigenfunction  solution  discussed  in  Referencps  2  through  4. 
These  are  compared  for  various  ka  values  as  well  as  incidence  angles.  The  limita¬ 
tions  of  the  PO  solution  are  also  discussed.  As  a  by-product  of  this  comparison, 
the  BSRCRCS  scattering  code  solution  for  the  disk  is  validated  against  an  exact 
eigenfunction  solution.  Another  issue  addressed  in  this  report  is  the  disk  curve 
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used  to  calibrate  a  bistatic  RCS  measurement  range.  The  disk  calibration  curve  is 
a  plot  of  the  specular  return  as  a  function  of  bistatic  angle.  This  curve  will  be  com¬ 
puted  for  both  horizontal  and  vertical  polarizations  and  various  ka  values  using  the 
PO  and  eigenfunction  solutions.  The  disk  has  become  a  possible  calibration  geometry 
for  bistatic  RCS  measurement  ranges  due  to  problems  encountered  using  the  sphere. 
The  bistatic  RCS  of  a  sphere  has  rapid  fluctuations  that  make  it  difficult  to  use 
as  a  calibration  standard  for  bistatic  radar  cross-section  measurements. 

The  derivation  of  the  PO  solution  is  given  in  Section  2.  The  PO  solution  cal¬ 
culates  the  scattered  field  using  the  geometrical  optics  (GO)  current,  where  the 
GO  current  is  denoted  by  Ufimtsev  as  the  uniform  current.  He  also  calculates 
the  field  scattered  due  to  the  non-uniform  current  on  the  edge  of  the  disk.  The 
superposition  of  fields  scattered  due  to  the  uniform  and  non-uniform  currents 
yields  the  total  scattered  field.  This  method  is  called  the  PTD.  The  BSRCRCS 
scattering  code  used  in  our  comparison  utilizes  PTD  as  its  method  of  solution; 
however,  it  is  not  a  closed  form  solution.  The  BSRCRCS  code  is  designed  fci 
arbitrarily  shaped  bodies  and  uses  numerical  techniques  to  perform  the  integrations 
required  to  implement  the  PTD  method. 

The  disk  is  a  coordinate  surface  in  the  oblate  spheroidal  coordinate  system;* 

therefore,  the  exact  eigenfunction  solution  can  be  found  using  the  degenerate  case 

of  the  oblate  spheroidal  functions  when  they  collapse  to  a  disk.  Detailed  discussions 

4 

of  this  approach  can  be  found  in  References  1,  2,  3,  and  4.  Mithouard  and  Hodge 
provide  a  computer  program  and  subroutines  required  to  calculate  the  RCS  based 
on  the  eigenfunction  solution.  This  computer  program  was  modified  by  Dominek  ‘ 
to  allow  its  use  for  higher  ka  values.  The  program,  as  modified  by  Dominek,  is 
used  in  this  analysis.  In  addition,  for  grazing  incidence  of  the  horizontal  polariza¬ 
tion,  this  computer  program  was  run  as  a  function  of  ka  to  determine  when  the 
eigenfunction  solution  is  no  longer  accurate.  This  was  done  to  provide  an  upper 
bound  on  our  use  of  the  eigenfunction  solution. 

The  details  of  the  PTD  and  eigenfunction  derivations  are  not  discussed  any 
further  in  this  report.  The  references  provided  give  detailed  analyses  of  these 
methods.  In  the  next  section,  we  derive  the  PO  based  expression  for  the  bistatic 
scattered  field  of  a  disk  illuminated  by  a  plane  wave. 


5.  Ufimtsev,  P.  la.  (1958)  Approximate  calculation  of  the  diffraction  of  a  plane 

electromagnetic  wave  by  certain  metal  objects,  II.  The  diffraction  bv  a 
disk  and  by  a  finite  cylinder,  Soviet  Physics  -  Technical  Physics, 

3(No.  1 1):2386-2396. 

6.  Dominek,  A.K.  personal  communication. 
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2.  BISTATIC  PHYSICAL  OPTICS  SOLUTION 


A  version  of  this  result  is  given  without  details  in  Reference  1,  and  the  details 
can  be  found  in  Reference  5.  The  following  result  was  derived  independently  be¬ 
cause  the  results  of  References  2  and  5  have  different  assumptions,  applicable 
geometry,  and  so  on,  than  the  eigenfunction  and  scattering  code  solutions;  there¬ 
fore,  it  was  decided  to  rederive  the  result  rather  than  to  make  the  requisite 
modifications  to  the  previous  results. 

The  plane  of  incidence  chosen  is  the  xz  plane  {ip1  =  0).  The  geometry  for  the 
disk  is  shown  in  Figure  1,  and  the  geometry  for  horizontal  and  vertical  polarizations 
are  shown  in  Figure  2  and  Figure  3.  For  our  given  plane  of  incidence,  the 
incident  E  and  H  fields  are  given  by 


Hl  =  Yq  r‘XEl  (lb) 

where 

Yo  =  l/Zo  =  sTmi  ,  (2a) 

k  =  277  /A  ,  (2b) 

k1  -  k  r1  =  -  k  (x  sin  4  +  z  cos  4 )  .  (2c) 

r1  =  -  x  sin  4  -  z  cos  4  >  (2d) 

9  1  =  x  cos  4  -  7  sin  4  (2e) 

01  =  Y  ,  (2f) 

and 

r  =  xx+yy+zz  .  (2g) 
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Figure  1.  Disk  Geometry 
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Figure  2.  Incident  Plane  Geometry  for 
Horizontal  Pol  (Pol  -  90)  |  |  Edge 


Using  this  current  (making  the  usual  far  field  approximations)  yields 
for  the  scattered  field  given  by 


ES  =  jk - 

po  r 


-j  k  r 


r  X  r  X  x 


r  X  r  X  y  cos 


D  (6>,  <p,  (,  ) 


whe  re 


a  277 

D  (0,0,5)=—  J  p'  f  k  [  sirl  0  c°s  <?>  -  V>’  >  h-  si”  ^  c°s 


0  0 


an  expression 


(5) 


'  dp'  .  (6) 
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P  (r,  0,o ) 


f  igure  .'I.  Incident  Plane  Geometry  for 
\  ertieal  I’ol  (l’ol  =  0)  Gdge 


We  cun  then  write  the  phase  term  of  the  integrand  in  the  form 
sin  8  cos  I o  -  <pM  +  sin  4  cos  9'  =  ft  cos  (•>  -9') 
such  that 


ft  cos  >  ^  sin  4  +  siii  f/  cos  0  . 


ft  sin  >  =  sin  8  sin  0  . 

■O' 


2  1 

ft.  =  'Gsin  8  sin  9  4  (sin  4  4  sin  6  cos  of” 


Therefore,  we  write 


a  2  it 

D(0.„.  $)  =  /  e^kp'ncos  ^  dp- 


(9) 


however 


2  77 

_L  I  eikp'ficosl  T  -f>  d,,,  =  J0(tlip.> 


(10) 


where  Jq  is  the  Bessel  function  of  order  zero.  In  addition. 


J.(kaft) 

J  p'  JQ(k ftp')  dp'  =  a  ~kT3 -  ‘ 


0 


(11) 


Hence,  the  PO  scattered  field  is  given  by 

■s  ;% 

■  9  cos  9  cos  0+0  sin  0 


_s  e*jkr  J  (kafl) 

ET  -  ja -  - 

PO  r  a 


cos  %  (0  cos  9  sin  0+0  cos  <p)) 


(12) 


We  want  to  identify  specific  input/output  polarization  pairs;  that  is,  what  is  the 
theta  polarized  scattered  field  for  a  phi  polarized  incident  field.  We  write  this  as 

g 

Therefore,  we  have  the  following; 


Eee  =  '  Ja“  r 


,-jk  <"  J,  (kafi) 


n 


cos  9  cos  q  , 


(13a) 


3  -  i  k  r  J.  (kafi) 

Kpe  s  J  a  r - 17 -  sin<p  • 


3  g-.i  k  r  (kaH) 

E  q  p  =  ‘  .i  a  — —  -  cos  4  cos  9  sin  tp  , 

s  g-i  k  r  J  ^  (k  a  O) 

E0«p  =  -  i a  r  — n — cos4  COS(^  * 


(13b) 


(13c) 


(13d) 
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These  expressions  can  be  modified  to  match  those  given  by  Ufimtsev'  for  the 
uniform  current. 


In  general,  the  scattering  cross  section  :s  defined  bv 


CT 


2 

lim  1  n  r“ 


(14) 


In  addition,  if  we  desire  the  polarization  dependent  cross  section. 
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we  define 


(15) 


where  a  and  b  represent  either  the  theta  or  phi  polarizations.  Therefore,  the 
input/output  polarization  dependent  cross  section  based  on  Eqs.  (15a)  to  (13d)  is 
given  by 


J  e  e  4 "  a 


J  t  (ka  S2) 

— n — 


2  2 
cos  6  cos  <j> 


a  =  4  7i  a ' 
00 


■  Ij  (k 

— n — 


.  2 

sin  p  , 


00 


4  a ' 


J  j  (ka  fi) 

- n - 


2  2  2 
cos  £  cos  0  sin  p 


'J  '4/7  a 
90 


.1  (kafi) 

n — 


2  2 
cos  £  cos”  p 


(16a) 


(16b) 


(16c) 


( 1 6  d ) 


It  should  be  noted  that  when  the  argument  of  the  Bessel  function  goes  to  zero  we 
have 


lim 


J  (ka«l 

— 52 - 


n  -  o 


k  a 

~T 


j0i°) 


k  a 

~T  ’ 


(17) 


There  are  some  special  cases  worth  mentioning. 

2  2  2 
a  =  «a  J  (2  ka  sin  0  )  cot  0  , 

and  for  the  specular  case, 

4  77  2  2 

<7  =  — j  (7  a  cos  0  ) 


In  the  case  of  backscatter, 

(18) 

(19) 
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This  has  the  proper  form  for  specular  scattering  given  by 


o 


spec 


(20) 


where  A  is  the  projected  area  of  the  scatterer.  Note,  both  the  backscatter  and 
specular  expressions  are  independent  of  polarization. 

These  expressions  were  programmed  and  used  for  comparison  with  both  the 
BSRCRCS  scattering  code  and  the  eigenfunction  code.  The  plots  in  Section  3  show 
where  the  bistatic  PO  solution  for  the  disk  breaks  down.  In  addition,  comparisons 
for  the  specular  response  using  the  PO  and  eigenfunction  solutions  will  be  shown. 
For  our  plane  of  incidence  (xz  plane),  horizontal  polarization  has  a  component 
parallel  to  the  edge;  therefore,  for  edge-on  incidence  (incidence  angle  nearly  90°), 
the  scattered  field  is  not  zero.  The  vertical  polarization  for  edge-on  incidence 
does  not  have  a  component  parallel  to  the  edge  and  hence  the  scattered  field  is  zero. 
The  eigenfunction  satisfies  the  proper  edge  condition;  therefore,  we  expect  the 
proper  behavior  at  the  edge.  This  will  be  evident  in  our  plots. 


3.  RESULTS 

In  this  section,  plots  that  compare  the  three  solutions  previously  described 
are  shown.  In  addition,  calibration  curves  computed  using  the  PO  and  eigen¬ 
function  solutiois  are  also  shown.  However,  before  beginning  the  comparisons 
it  is  of  interest  to  determine  when  the  eigenfunction  solution  breaks  down  as  a 
function  of  the  electrical  circumference  of  the  disk  (ka).  For  horizontal  polariza¬ 
tion,  we  plotted  backscatter  and  forward  scatter  for  grazing  incidence.  As 
expected,  the  backscatter  case  accounts  for  the  creeping  wave  around  the  disk 
edge;  hence,  we  see  the  interference  lobes.  On  the  other  hand,  the  forward 
scatter  is  not  affected  by  this  mechanism  and  is  therefore  a  smoother  curve. 
Figure  4  shows  this  behavior.  Also,  we  can  see  that  the  solution  breaks  down 
beyond  ka  =  22.  Therefore,  we  only  used  this  solution  up  to  ka  -  20. 
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SCATTERING  FROM  A  DISK  (CRAZING  INC:  HORIZ  POL) 


Figure  5  through  Figure  B  show  the  bistatic  radar  cross-section  (RCS)  of  a 
disk  illuminated  by  a  plane  wave  at  normal  incidence.  The  plots  shown  are  for 
both  horizontal  and  vertical  polarization  at  several  ka  values.  For  all  these  cases 
the  phi  observation  angle  is  zero  (same  as  phi  incident).  For  horizontal  polariza¬ 
tion  (Pol  =  90),  the  agreement  is  better  than  for  vertical  polarization  (Pol  =  0). 

As  the  disk  ka  increases  and  the  observation  angle  exceeds  45°,  the  PO  solution 
begins  to  break  down;  however,  the  agreement  between  the  BSRCRCS  code  and 
the  eigenfunction  code  remains  reasonable  for  larger  angles,  especially  for 
horizontal  polarization. 
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Figure  5.  Bistatic  Scattering:  ka  =  5:  Theta  Inc  =  0: 
Horiz  Pol:  Phi  Inc  =  Phi  Scat  =  0 
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Figure  6.  Bistatic  Scattering:  ka  =  5:  Theta  Inc  = 
Vert  Pol:  Phi  Inc  =  Phi  Scat  =  0 
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Figure  8.  Bistatic  Scattering:  ka  =  10:  Theta  Inc  =  0: 
Vert  Pol:  Phi  Inc  =  Phi  Scat  =  0 
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Letting  p  =  180°,  we  are  in  the  quadrant  where  the  specular  lies  for  our  given 
plane  of  incidence.  We  see  that  regardless  of  how  different  the  results  are  in 
general,  they  agree  (within  a  few  dB)  at  specular.  Figures  9  through  14  show 
these  results  for  several  incident  angles.  The  agreement  at  specular  gives  us 
hope  the  calibration  curve  will  agree  over  a  wide  range  of  bistatic  angles. 

Plotting  the  calibration  curve,  we  see  that  the  calibration  curve  calculated 
using  PO  goes  to  zero  for  both  polarizations;  however,  the  calibration  curve 
calculated  using  the  eigenfunction,  as  expected,  does  not.  We  see  reasonable  agree¬ 
ment  for  both  polarizations  over  bistatic  angles  less  than  140  degrees.  Figure  15 
through  Figure  20  shows  this  comparison  for  several  ka  values.  As  ka  increases, 
the  agreement  between  the  PO  and  eigenfunction  solutions  for  specular  improves. 

As  we  see  in  Figure  20,  ka  =  20  the  difference  is  negligible  for  0  £  |3  S  140°. 

4.  CONCLUSIONS 


Our  eigenfunction  code  is  only  accurate  up  to  ka  =  20.  Beyond  that  it  is  no 
longer  useful. 

We  saw  that  the  PO  solution  for  normal  incidence  agreed  very  well  with  the 
eigenfunction  and  BSRCRCS  codes  over  the  main  lobe  and  a  significant  number  of 
side  lobes.  For  small  ka  the  agreement  is  remarkable. 

When  we  concentrate  on  the  specular  return  from  the  disk,  PO  docs  a  good 
job  for  incidence  angles  as  large  as  70  degrees.  This  is  really  not  that  surprising 
since  the  specular  return  is  an  optics  quantity  governed  by  geometrical  optics.  In 
fact,  the  stationary  phase  evaluation  of  the  PO  integral  gives,  in  general,  the  GO 
reflected  field.  As  we  saw  in  the  calibration  curves  (specular  response),  the  PO 
and  eigenfunction  solutions  agree  very  well  for  large  ka;  however,  for  small  ka 
we  would  need  a  correction  factor  to  use  the  PO  curve  for  calibration.  We  saw 
that  the  disk  calibration  curves  are  well  behaved  and  have  no  rapid  fluctuations  and 
there  are  several  available  methods  to  calculate  the  curve.  This  makes  the  disk 
an  ideal  candidate  as  a  standard  for  bistatic  RCS  measurements.  We  also  saw  that 
the  PO  and  eigenfunction  solutions  generate  the  same  curve  over  a  wide  range  of 
bistatic  angles  (0  <  (3  <  140°)  for  large  ka.  This  means  there  exists  an  easy  method 
of  generating  the  applicable  calibration  curve  for  a  particular  disk  as  long  as  one 
remains  within  the  valid  angular  region  with  an  appropriately  sized  disk. 


13 


Figure  9.  Bistatic  Scattering:  ka  =  10:  Theta  Inc  =  15: 
Horiz  Pol:  Phi  Inc  =  0:  Phi  Scat  =  180 
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Figure  10.  Bistatic  Scattering:  ka  =  10:  Theta  Inc  =  15 
Vert  Pol:  Phi  Inc  =  0:  Phi  Scat  =  180 
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Figure  11.  Bistatic  Scattering:  ka  =  10:  Theta  Inc  =  45: 
Horiz  Pol:  Phi  Inc  =  0:  Phi  Scat  =  180 
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Figure  12.  Bistatic  Scattering:  ka  =  10:  Theta  Inc  =  45: 
Vert  Pol:  Phi  Inc  =  0:  Phi  Scat  =  180 


5 


O  CD  10  El 


SPECULAR  SCATTERING  FROM  A  DISK  (CO-POL) 


Figure  13.  Specular  Scattering:  ka  =  5:  Horiz  Pol 
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Figure  16.  Specular  Scattering:  ka  =  5:  Vert  Pol 
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